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We present preliminary results for the decay constant of the Bs meson, /s^, at three values of f3 — 5.7, 6.0 and 
6.2 using NRQCD and clover fermions for the heavy and light quarks respectively. As a consistency check the 
decay constant has also been extracted from the axial-vector current at finite momentum. In addition, we discuss 
the cancellation of 0{a/ (aAIo)) terms and the remaining uncertainty in /b, from higher order divergences. 



1. Simulation details 

For the heavy quark we use an NRQCD action 
consistent to 0(1/Mg): 

QHAt + Ho 
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where, 
Ho = 

SH = 
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The 0(1/Mq) correction to the kinetic en- 
ergy (expected to be the largest contribution from 
this order) and the first two discretisation cor- 
rections are also included. We implement tad- 
pole improvement throughout, using the plaque- 
tte definition of uq, and set c; = 1. 

For the light quark we use the clover action 
with the tadpole improved value of csw at /3 = 
5.7 and 6.0 and the non-perturbatively deter- 
mined value at /3 — 6.2; csw ~ 1-61 as deter- 
mined by the Alpha collaboration [Q, compared 
to Csw — 1-48 using tadpole- improvement. The 
configurations at /3 = 5.7 and the configurations 
and light quark propagators at f3 — 6.2 were gen- 
erously provided by the UKQCD collaboration. 

*In collaboration with C. T. H. Davies and J. Hein, Glas- 
gow Univ.; J. Shigemitsu, Ohio-state Univ.; C. Morn- 
ingstar, San Diego; A. Ali Khan, Tsukuba. 



f3 


V 


N 


a-' {nip) (GeV) 




5.7 


12^ X 24 


278 


1.14(3) 


- 4.2 


6.0 


16^ X 48 


102 


1.92(7) 


- 2.2 


6.2 


24^ X 48 


144 


2.63(9) 


~ 1.6 



Table 1 

The simulation details. The errors on a^^ include 
statistical errors and those due to the chiral ex- 
trapolation of Trip. N denotes the number of con- 
figurations. 



The light quark mass is fixed to the strange quark 
mass using the K meson mass, with the uncer- 
tainty in this determination estimated by fixing 
TTiq using the 0. Further simulation details for the 
3 ensembles are given in table 0. 

The pseudoscalar decay constant is defined as 



\Ap\PS) 



QCD 



(4) 



in Euclidean space. On the lattice matching fac- 
tors Ci relate the lattice operators to the current 
in full QCD. For the zeroth component of the cur- 
rent to 0{a/M): 



{Aq)qcd = 



(5) 



where. 
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Figure 1. Preliminary results for /b^ as a func- 
tion of the lattice spacing. Our results (circles) 
are compared to those from JLQCD (squares) ||] . 
All errors include statistical and systematic un- 
certainties added in quadrature. 

The 0{aa) discretisation error in the current is 
removed by defining 1^, 

j0.^mp ^ jf^^ +CAjf""^ (9) 

jf-) = aq{rhh5loQ (10) 

The CiS have been calculated to 1-loop how- 
ever, the q* (which also depends on aMp) at 
which the strong coupling, a, is computed is not 
yet known. Thus, we average the results obtained 
using aq* = 1.0 and tt. Note that in the static 
limit aq* ~ 2 ||^ for Wilson light fermions. 

The results at /3 = 5.7 and 6.0 have appeared 
previously in |^ and |^ respectively. Further de- 
tails of our methods and analysis can be found in 
these references. 

2. Scaling of fs. 

Our values for fs , calculated to 0{a/M), are 
presented in figure |l| as a function of a~^. The 
results are consistent with scaling, and /b^ ~ 
200 MeV. Nice agreement is found with the re- 
sults of the JLQCD collaboration ||^, also shown 
in the figure. This group uses the string tension 
to set a, and an O(l/Mo) NRQCD action. How- 
ever, this is unlikely to affect the comparison sig- 
nificantly. 

Table details how the errors are estimated. 



Source 


13 = 5.7 


6.0 


6.2 


statistical 


3 


3 


2 


disc. OiiaAQCDf) 
pert. 0(a2,a/(aM)2) 
NRQCD 0(1/M2) 


13 


4 


2 


13 


8 


9 


1 


1 


1 


Kg 


+4 


+4 


+4 




3 


4 


3 


Total 


19 


11 


11 



Table 2 

Estimates of the statistical and main systematic 
errors, in percent, in our values for Jb^ 



Note that - 1-2.6 GeV is the range in which 
the NRQCD action and clover light fermions can 
be applied to the B meson. For coarser values of 
a, the discretisation errors from the light quark 
action (and gauge action) increase rapidly, as do 
the O(a^) perturbative errors. Conversely, if a 
becomes too fine, the discretisation errors are un- 
der control but the 0(a/(aM)^) perturbative er- 
rors increase dramatically as aMg drops below 
1. Overall, the perturbative errors are the main 
source of uncertainty, although at /3 = 5.7 the 
discretisation errors are of equal magnitude. 

3. Jb^ extracted at finite momentum 

In order to investigate momentum dependent 
discretisation errors in the decay constant we 
computed the ratio 

{jf\^/{jf\=^E{p)jM (11) 

where {j]j^^)p = f'^^jE{'p) (without rcnormalisa- 
tion). The RHS of equation |ll| is a slowly varying 
function of \p\ , which is close to 1 for the range of 
momenta we studied (up to 1.5 GeV ai (3 — 6.2). 
Our results, presented in table are in agree- 
ment with this expectation at /3 = 6.0 and 6.2, i.e. 
the momentum dependent discretisation errors in 
Jb^ are not significant. However, at /3 = 5.7 a 
10—20% deviation from 1, is seen as \p\ increases, 
although this is within the magnitude expected 
for 0((ap)^) discretisation errors. 

4. Power-law divergences 

Matrix elements in NRQCD beyond zeroth or- 
der diverge as aA/o 0. In the case of /b 
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{fwm)/{fWM) 


a 


= 1 


2 


3 


4 


5.7 


0.94(1) 


0.88(1) 


0.84(1) 


0.82(2) 


6.0 


1.01(1) 


1.01(3) 






6.2 


1.00(2) 


1.02(3) 


1.03(4) 


1.07(6) 



Table 3 

The decay constant extracted at finite momentum 
for aAfo close to aM^. |p| = 2mT/{aL), where 
n = 0,1, . . . and L is the spatial extent of the 
lattice. 

there are unphysical, ultra-violet, contributions 
to (J|,'^),i > 0, which are cancelled order by order 
in perturbation theory by terms appearing in the 
perturbative coefficients. In general, simulations 
are performed using aAfo>l and so the unphysi- 
cal contributions are not expect to be large and 
their cancellation should be under control. 

Considering equation |^ in more detail (see 
for definitions). 



{Ao)qcd = (l + apo)(jf ) + 



+apdisc{Jl ) 

The explicit contributions to po are 
Po = [-Bo - 2 (C? + C'q) ~ Coo - Cio]- 



(12) 



(13) 



The lowest order divergent contribution to the 
current is 0{a/ {aM)), which appears through 
the tree- level term {j')^'^) and is cancelled by 
the mixing term aC,io{jf'^) , where ^iq is the 
renormalisation due to the mixing between J^^^ 

and jj,^"*. The remaining divergent contributions, 
0{a/{aMY,a'^ /[aM)) etc, appearing in equa- 
tion are cancelled at higher orders; the un- 
certainty in /s^ due to these remaining terms is 
well within our estimates of the systematic uncer- 
tainties in table ||. 

In fact a large part of (j},^'') is unphysical, 

aCio(4" V(4^^> ^ 0-5 - 0.6 for aU /3s and aMo 
for aq* = 2.0. However, once the 0{a/{aM)) 

contribution to (j'}^^) is cancelled, the scaling 
behaviour of this term improves, as shown in 
figure |2| suggesting the remainder is physical. 
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Figure 2. The contribution to the decay constant 
from the tree-level current correction before and 

),ag* = 2.0. 



after the subraction of the q;Cio(«/1 



aCio(4°^) - 0.025GeV3/2 is 



However, (J^^^) 

small, (cf fVid ~ 0.500GeV^/^), and of the or- 
der of the higher order perturbative and NRQCD 
corrections. Hence, we cannot determine the size 
of the physical part of {J^'') reliably. We empha- 
sise this does not lead to a significant uncertainty 
in /b^ , nor in the slope of the decay constant with 
1/Af, which are dominated by (J^^^) and the cor- 
responding perturbative coefficient (without the 
Cio term). 
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